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Abstract 

Shape-invariant signals under Fourier transform are investigated leading to a class of eigenfunctions for the Fourier operator. The 
classical uncertainty Gabor-Heisenberg principle is revisited and the concept of isoresolution in joint time-frequency analysis is intro¬ 
duced. It is shown that any Fourier eigenfunction achieve isoresolution. It is shown that an isoresolution wavelet can be derived from 
each known wavelet family by a suitable scaling. 
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1 Preliminaries 

The Fourier transform is often interpreted as a linear operator &. 
An interesting problem in this framework is to find out the eigen¬ 
functions in the language of operators CHS- Let V be a vector 
space equipped with a linear transform, T : Y —> Y, v H > T(v). 
Under the linear transform T, eigenfunctions are solutions of 
T(\) = X - which corresponds here to J F{f(t)}(co) = X •/(©) 
where / € L 2 (K.) and X is a scalar. They are a quite remark¬ 
able class of functions, which preserves the shape under Fourier 
transform: Both the signal and its spectrum (time and frequency 
representation) have the same shape. In joint time-frequency rep¬ 
resentation 00 this feature can represent a very good balance 
between the two domains. It is well known that the Gaussian 
pulse is a signal whose shape is preserved under the Fourier op¬ 
erator: 

e -d/2^y^. e -® 2 /2. 

This can easily be derived by writing 

-L e - ,2 / 2 -e J0} 'dt =F{(0). 

V2n J— 

Deriving this equation and using integral by parts, one notice 
that: = —a>F((o). The solution of the differential equa- 

tion ^F(co) + COF(co) = 0 under the initial condition F(0) = 1 
is F (to) = e~“ 2 / 2 . It follows promptly that X = \[2%. 
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The question is: Are there other eigenfunctions? This mat¬ 
ter is addressed in the next section. It is worthwhile to bear in 
mind that some results in this paper are deliberately non nova, 
sed nove. 

2 Shape-invariant Signals: Eigenfunctions of the 
Fourier Operator 

Let £{•} and 0{} denote the functionals that extract the even 
and odd part of a given signal, respectively. 

Proposition 1 Let /(f) < — >F(co) be an arbitrary Fourier trans¬ 
form pair. Then the signal 

h{t) = V2 k • £{/(f)} + £{F(®)} 

is invariant under the Fourier transform. Furthermore, we have 
that: H(co) =&{h{t)}{(o) = s/2n-h{(0). 

Proof: It follows from the definition of h(-) that 

2 ■ h{t) = V2 n-\f{t)+fi-t)\ + [F{t) +F(-t)}. 

Taking the Fourier transform, 

2-Hico) = V2k- [F((o) +F(-g>)] + [2nfi-(o) + 2nf{(o)\. 

and the proof follows. □ 

Corollary 1 Each even function f(t) <—> F(o)) induces a 
Fourier invariant h(t) = \/2nf(t) + F(t). 

For instance, the following signals 

hi(t) = V2n • - 1 , + ne ^, 

1 +t 

,— 2 

/ 12 (f) = v2n\t\ - — 

t L 
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have spectra with similar shape. Another remarkable example is: 



where sech(-) is the hyperbolic secant function. 

Proposition 2 Let f(t) < — >F(co) be an arbitrary Fourier trans¬ 
form pair. Then the signal 


Thus, — co 2 F(co) = kXF{(o), i.e., its spectrum also 

obeys a similar differential equation. Therefore, / and F have 
identical shape, since they are solutions of the same differential 
equation. □ 

The key equation for shape-invariant signal is thus 4^f(t) — 
t 2 f(t ) = Kf(t '). Let us try solutions of the form 

f{t)=p{t)e~ t2/2 , 


h(t) = V2ii-Q {f{t)} - 0 {F(f)} 


where p(t) is a function to be determined. Therefore, 


is an invariant under Fourier transform. Furthermore, 

fiF {/;(/)} = —\/2nh{(d). 

Proof: The proof is similar to the proof of Proposition [I] 

Corollary 2 Each odd function f(t ) <—> F(co) induces a 
Fourier invariant h(t) = \/2nf[t) — F(t). 

Let us now focus on a particular and important class of Fourier 
invariant, which generates an orthogonal and complete set. To 
begin with, let us denote by S the class of eigenfunctions of the 
Fourier operator defined according to the following proposition. 


At 2 f 


P(t) 


,~t 2 / 2 


— t 2 p(t)e ?2 / 2 = Kp(t)e ' 2 / 2 . 


After simple algebraic manipulations, we derive 
d 2 d 

di 2 ^ ~ 2t di P ^ + ( K+1 ' )p ( t ' ) = °’ 

where n is a integer. 

A standard differential equation of the above form (§1 is 

d 2 d 

- 2t—p(t) +2np{t) = 0 , 


( 2 ) 


Proposition 3 A signal f(t) is in $ if and only if, the signal f 
satisfies the differential equation 4-zf(t) — t 2 ■ f(t ) = K ■ f(t), for 
some scalar KG C. 

Proof: We begin demonstrating the suffiency. By hypothesis, we 
have that 

f(t)^Xf{co). 

The properties of time and frequency differentiation for -F 
give: 

^ (7®) 2 A/(0>), 
dta- 

Adding above expressionfQ. we derive 

^/(®)-® 2 /(co) . 

Thus, the signal ^f(t) — t 2 f{t ) has also its shape preserved, 
provided that / itself preserves its shape. Therefore, j^rfit) — 
t 2 f{t) G S, that is, we are looking for signals such that ^3 f(t) — 
t 2 f{t) = Kf{t), since they have identical eigenvalues. 

Now we demonstrate the necessity. 

By hypohthesis, the signal fit) satisfies the differential equa¬ 
tion jp;f(t) — t 2 f[t) = Kf(t), K G C. Applying the operator ,F, 
we obtain: 

7 d 2 

Uco) 2 F{(o) + ^2^0) = kAF(o). 

'N.B. Subtracting: f( t )+ t 2 f(t) FF -X [-£rf{(o) + ffl 2 /(to)]. 


where n is a integer. Thus, for a suitable choice K = — (2 n + 1) 
(eigenvalues), the solutions p(t) are exactly Hermite polynomi¬ 
als ©, which form a complete orthogonal system. Thus, we 
have: 

p(t)=H„(t), 

where 

H 0 (t) = \, 

H\ (?) = 2t, 

Hi{t) = —2 + 4r, 

// 3 (f) = -12f + 8f 3 , 

Hfit) = 12 —48f 2 + 16f 4 , 


Proposition 4 Possible eigenvalues of the Fourier transform are 
the four roots of the unit (± 1 , ± 7 ) times y/2n. 

Proof: Let us denote by the operator corresponding to it- 

_ jr ^ yp 

erate n times the operator ,F. Let t <—> co <—> co <—7 Q. be 
the Fourier domain variables for the iterate Fourier transform. 
Observe that, for f G S’, we have: 

JF^{f(t)}(Ll)=4n 2 m). 

But, 

^ ( 2 ) {/(0}(a)') = A 2 /(-co'), 

&W{f(t)}(Q.) = X\m. 

From ([3]) and ©, it follows that A /y/27C G C has order 4. □ 


A -/ N t 2 r ,.s .&. 
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s 2 / 'i 00 

We conclude that < \if n (t) = H n (t)e~ r ~' 2 > are shape- 
l J n =0 

invariant under Fourier operator associated to A„ = (—j) n y/2n. 
Therefore, 

H n (t)e- ,2 S 2 A- (-j) n V2tH n (co)e-“> 2 / 2 . (5) 

Another interpretation can be derived evoking Rodrigues’ for¬ 
mula l 6 l : 

H n (t) = (-\) n e' 2 ^e~' 2 . 

The 2nd-order differential equation hold by invariant signals is 

d 2 2 \ 

-^y + (2 n + 1 - x 2 )y = 0. 

The above differential equation is exactly the celebrated 
Schrodinger equation for the harmonic oscillator Q. 

3 Consequences on the Time-Frequency Plane 

Let us now investigate certain consequences of eigenfunctions of 
the Fourier operator on the time-frequency plane mm. 

Let f(t ) be a finite energy signal E, equipped with Fourier 
transform, F(co). The time and frequency moments of / are de¬ 
fined by: 

pr r^ritrmdt 

nrw/(odr ’ 

= yy\m\ 2 « t 

w _ r <x ,F*(co)co n F{co)dco 
^ f-ooF*(a>)F((0) dto 

By analogy to Probability Theory, the term \f(t)\ 2 /E denotes 
a “time-domain” energy density, where £ is a normalising factor 
so as to make the whole integral of the density be equal to the 
unity. It is customary to deal with the energy spectral density 
G(co) = \F(o))\ 2 , whose integral over a frequency band gives the 
energy content of the signal within such a band. Let us suppose 
in the sequel, without loss of generality, that E = 1 (energy nor¬ 
malised signals). 

The “effective duration” (respectively the “effective frequency 
width”) of a signal /(f) (respectively F(co)) is defined according 
to: 

Ar = <\J2n(t — t) 2 r.m.s duration, 

A/ = \/2n(f — f) 2 r.m.s bandwidth, 

where Ar and A/ correspond to the standard deviation, i.e., 
spreading measures. However, other common and much hand¬ 
ier definitions are 

A r = \/V- ? ) 2 ’ 

A ffl = 

Clearly, A, = At/y/2K and A ffl = \f2%Af. 


3.1 Revisiting the Gabor Principle 

By applying arguments from quantum mechanics (7), Gabor j9] 
flOl derived an uncertainty relation nowadays called Gabor- 
Heisenberg principle for signals: A t A/ > 1/2, proving that 
time and frequency cannot be exactly measured (simultane¬ 
ously). The Gabor-Heisenberg uncertainty principle states a 
lower bound on the product Ar ■ Act), or alternatively: 

A,-A* >1/2. (6) 


Proposition 5 The Gabor lower bound is only achieved by the 
first invariant signal (eigen functions of J? operator). 


Sketch of the proof: From ([6]), the bound is achieved if, and only 
if, -tj f(t) = tctf(t) This condition can be interpreted as: “deriva¬ 
tive in time domain” is equivalent to the “’derivative in frequency 
domain”. Therefore, 


57 T/W = * 




Kf(t) + K 2 t 2 f(t). 


Simple manipulations yield: 

^2 /(0 ” K ’( 1 + Kf 2 )^) 2 -fit) = 0 . 

The only solutions on $ correspond to K = ±1, i.e., 4^f(t) + 
(1 - r 2 )/(r) = 0 or £/(r) - (1 +r 2 )/(r) = 0. 


Proposition 6 Any real signal f(t) <—> £(m) such that 
f(t),jjf(t),F(co), ^F(co) e L 2 (R) have finite resolutions. 

Proof: Applying the Parseval-Plancherel Theorem (6), it follows 
that 

J t 2 f 2 {t)dt = J [jtf(t)} ■ [jtf{t)]*dt 


ir 

2k J c 


dm 


F(co) 


and 


J co 2 |£ 2 (m)|“dm = J [jcoF((o)] ■ \jcoF((o)Y 


dm 


= 2n 


j: 


12 


d t 


m 


Therefore, 


Af = 


/-~ld^(®)l 2d ® 


f -00 |£(®)| 2 dm 


2 _ f-00 1 37/(0 |~ df 

m/(oi 2 dco 


A^) = 


< 00 , 


< OO. 


Thus, the above quantities are given by the square root of the 
ratio between the energy of the signal derivative and the energy 
of signal itself. Thus, the resolution for the Fourier invariant 
signal sech(-) given by (U} is 


Af — A® — 



0.7235987766... 
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since that 


J sech(/) 


/; 
/:.(!) 


dt = 2, 
2 { 


tanfr(/)sech"(f)dr = 

2 


seclr(f)dr = 


where tanh(-) is the hyperbolic tangent function. 


The quantity y/ A t / A<d is refered to as the isoresolution factor. 
The essential idea of isoresolution can be placed in the wavelet 
structure. Normally, the basic wavelet of a family — 

vM 

holds the admissibility condition but often does not achieve 
isoresolution. We propose here to redefine the basic wavelet of a 
family so as to achieve isoresolution. For instance, the standard 
Mexican hat wavelet i/Mhat^) satisfies: 



2 (r 


l)' 


/ 2 

y/nV 3 


& 



Proposition 7 ([9j) Time-frequency uncertainty of Fourier 
Eigenfunctions \ff *(?) = H n (f)e~ t2 ! 2 e^ 1Jr ^, where (Oq and (/)o 
are constants, attain quantized values of the Gabor-Heisenberg 
lower bound, i.e. 

At • A/ = i • (2n + 1), 

Ar ■ A ffl = 2 ' (2w + 1). 

That is why Gabor functions are relevant in some problems 

(e.g. HD.) 

4 The Concept of Isoresolution Wavelet 

The concept of isoresolution analysis is introduced in this sec¬ 
tion. According to the Gabor principle, if one increases reso¬ 
lution in one domain, the resolution must decrease in the other 
domain so as to guarantee the lower bound given by ©. When 
analysing signals in joint time-frequency plane, frequently, there 
is no grounds to assure a better resolution in a domain than in 
the other domain. As an interesting property, any Fourier eigen¬ 
function achieves isoresolution as it can be seen by the following 
proposition. 

Proposition 8 Fourier-invariant signals perform an isoresolu¬ 
tion, that is, At = A ffl 

Proof: Supposing that f £ S, then F(tx>) = X f(co). Therefore: 

ff' 00 F(co)co 2 F(co)* dco ff 00 co 2 \X\ 2 f 2 (co)dco 
jrjF(co)\ 2 dco = JfJX\ 2 f 2 (co)dco 


The isoresolution Mexican hat wavelet can be found applying 
Proposition [9j 



For any isoresolution wavelet, the scaling by a > 1 or a < 1 cor¬ 
responds to unbalance resolution in a different way. Table Q]dis- 
plays both time and frequency resolution for a few known contin¬ 
uous wavelets: Gaussian derivatives, Mexican hat, Morlet, fre¬ 
quency B-Spline, Shannon, and Haar IfTZl . The wavelet Gausl 
is an invariant wavelet therefore it achieves isoresolution, in ac¬ 
cordance to proposition 8. It is valuable to mention that compact 
support wavelets (in time or frequency) cannot attain isoresolu¬ 
tion, since no signal can simultaneously be time and frequency 
limited |13l . 

5 Perspectives and Closing Remarks 

Eigenfunctions of the Fourier operator were investigated and the 
Gabor principle was revisited defining the concept of isoreso¬ 
lution, i.e, a signal with the same time and frequency resolu¬ 
tion. The functions {y/,,(f)} ( see ©) turn up as a very ap¬ 
pealing choice for designing representations such as wavelets. 
It is time to try finding new wavelets starting with (0. Since 
they are solutions of a wave equation (2nd order differential 
equation), our approach (Mathieu ED, Legendre ED, Cheby- 
shev |T6)) can be useful to construct new wavelets: The Quan¬ 
tum Wavelets, or Gabor-Schrodinger wavelets. The construction 
of new wavelets based on these complete, orthogonal, domain 
shape-invariant system is currently being investigated. The idea 
is to adapt the concept of isoresolution in orthogonal multireso¬ 
lution analysis Il7lll8l . 


and the proof follows. □ 

This is an interesting property for signalling on the joint time- 
frequency plane. 

It is suggested here the changing of the time-frequency reso¬ 
lution by a proper scaling that allows for identical resolution in 
both domains. 
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Proposition 9 If y/(t) has effective duration At and effective 
bandwidth Aft), then the scaled version t// (^y/ A, / A m t^j achieves 
isoresolution. 

Proof: Scaled versions y/(at), a f 0, have resolutions A,/|a| and 
|a| • A ffl , so |a| can be appropriately chosen. □ 
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